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THE SPACE DENSITY OF COMPTON THICK AGN AND THE X-RAY BACKGROUND 

E. Treister^-^ '*, C. Megan Urry^ "''*' and Shanil Virani''''' 
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ABSTRACT 

We constrain the number density and evolution of Compton-thick Active Galactic Nuclei (AGN) . In 
the local Universe we use the wide area surveys from the Swift and INTEGRAL satellites, while for liigh 
redshifts we explore candidate selections based on a combination of X-ray and mid-IR parameters. We 
find a significantly lower space density of Compton-thick AGN in the local Universe than expected 
from published AGN population synthesis models to explain the X-ray background. This can be 
explained by the numerous degeneracies in the parameters of those models; we use the high-energy 
surveys described here to remove those degeneracies. We show that only direct observations of CT 
AGN can currently constrain the number of heavily-obscured supermassive black holes. At high 
redshift, the inclusion of IR-selected Compton-thick AGN candidates leads to a much higher space 
density, implying (a) a different (steeper) evolution for these sources compared to less-obscured AGN, 
(b) that the IR selection includes a large number of interlopers, and/or (c) that there is a large number 
of reflection-dominated AGN missed in the INTEGRAL and Swift observations. The contribution of 
CT AGN to the X-ray background is small, '^9%, with a comparable contribution to the total cosmic 
accretion, unless reflection-dominated CT AGN significantly outnumber transmission-dominated CT 
AGN, in which case their contribution can be much higher. Using estimates derived here for the 
accretion luminosity over cosmic time we estimate the local mass density in supermassive black holes 
and find a good agreement with available constraints for an accretion efficiency of ~10%. Transmission- 
dominated CT AGN contribute only ~8% to total black hole growth. 
Subject headings: galaxies: active, Seyfert; X-rays: galaxies, diffuse background 



1. INTRODUCTION 

It is now clear that most accretion of mass onto 
super massive black holes is obscured from our view 
(e.g. , iFabianI 1X9991 : iTreister et all l2004f ) . Observations 
of nearby Active Galactic Nuclei (AGN) suggested that 
the lo cal ratio of obscured to unobscured sources is ^4:1 
(e.g. iRisaliti et al.llT999f l. A similarly high fraction of 
obscured AGN has been used to explain the spectrum 
and normalization of the extragalactic X-ray Background 
(XRB), as shown by the latest AGN populat ion synthe- 
sis models (|Treister fc 1177^120051 : iGilli et al.l[2007 ). The 
XRB gives an integral constraint to the AGN population 
and its evolution; the most recent deep surveys show that 
^90% of the observed 2-8 k eV XRB radiation can be at- 
tributed to resolved AGN (|Hickox fc MarkevitchI l2006l 
and references therein). 

The most obscured AGN known are those in which 
the neutral hydrogen column density {Nh) in the line of 
sight is higher than the inverse Thomson cross section, 
A''//~1.5x 10^^ cm~^. These are the so-called Compton- 
thick (CT) AGN. If the obscuring column density is 
smaller than ~10^^ cm~^, direct emission from the nu- 
cleus is still visible at energies greater than ~10 keV, 
while the radiation at lower energies is completely ob- 
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scured by photoelectric absorption; in this case we have 
a transmission-dominated CT AGN. For sources with 
A''^f>10^^ cm~^ the X-ray emission is significantly af- 
fected by Compton scattering at all energies, fully ob- 
scuring the direct AGN emission and leaving only the 
much fainter reflection component to be detected; these 
are reflection-dominated AGN. 

Contrary to the situation for less obscured sources, 
not much is known about the number density of CT 
AGN. Thanks to the deep Chandra and XMM-Newton 
surveys it is now clear that the fraction of moderately 
obscured, Compton-thin, AGN is on average ^3/4 of all 
AGN, and is higher a t lowe r luminosities (Ueda et al.l 
120031 ITreister fc Uml 120051: iBa rgcr et al. 2005 ) and 
higher redshifts (iLa Franca et al.ll2005i : .Ballantyne et al.l 
120061 ITreister fc Urrvll2006t l. but there are no compa- 
rable constraints on the number of CT AGN. About a 
dozen CT AGN have been ide ntified in the local Uni- 
verse ([Delia Ceca et all l2008al and references therein) . 
In fact, two of the three nearest AGN ar e Compton 
thick (NGC4945 and the Circinus Galaxy; iMatt et all 
I2OOOI). Based on a s ample of 49 local Seyfert 2 galax- 
ies, iGuainazzi et al.l (|2005f ) estimated that ~50% of all 
obscured AGN (7Vjy>10 cm~^) are Comptoii thick , 
and similar estimates were made bv IRisaliti etall (|1999f ) 
based on much smaller numbers. 

However, so far there has been no systematic study of 
the statistical properties of CT AGN with a well-defined 
selection function. Hence, while it has been hypothesized 
that CT AGN ca n contribute up to ~30% of the XRB 
(|Gilli et al.l |2007( ) and represent a significant fraction 
of the cosmic accret ion onto supermassive black holes 
(Mar coni et al.l 120041 ). this has not been demonstrated. 
Now, thanks to the wide-area surveys at high energies 
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perfo r med with the INTEG RAL /IBIS ()Beckmann et all 
120061 iKrivonos etlll I2007D an d the Swift/Burst Alert 
Telescope fBAT: iTueller et al.ll2008[ ). it is possible to 
study a well-defined sample of CT AGN in the local Uni- 
verse. Furthermore, since most of the absorbed energy 
is re-emitted at mid-IR wavelengths, deep observations 
with the Spitzer observatory can be used to select CT 
AGN candidates at high redshift, z~2 (iDaddi et al.ll2007t 
iFiore et al.ll2008t I Alexander et al.ll2008D . vieldii^n up- 
per limit to the number density of these sources. 

In this paper we constrain the number density of CT 
AGN in the local Universe from high-energy observa- 
tions, and at high redshift using a combination of X-ray 
and mid-IR data. We compare the observed numbers of 
CT AGN with expectations from AGN population syn- 
thesis models that explain the XRB emission and we 
study the degeneracies affecting these models. Finally, 
we compute the implied density of supermassive black 
holes as a function of redshift, including transmission- 
dominated CT accretion. When required, we assume a 
ACDM cosmology with /io=0.7, and f]A=0.7, 

in agr eement with the mo st recent cosmological observa- 
tions (jSpergel et al.l 120071 ). 

2. THE LOCAL SAMPLE OF CT AGN 

One of the best ways to find CT AGN is by observ- 
ing at high energies, namely E>10 keV. The hard X- 
ray spectrum of a CT AGN is characterized by at least 
three components: an absorb ed power law with an up - 
per cutoff at ~300 keV (e.g., iNandra fc Poundsl 11994 ). 
a Compton reflection hump which peaks at ~30 keV 
(|Magdziarz fc Zdziarskil |1995[ ) and an iron Ka line at 
~6.4 keV. Not all components are clearly observed i n 
all AGN fe.g.. ISoldi et~an [20051 : IFeckmann et al.ll2004l ). 
perhaps because of the low signal-to-noise of some of the 
observations. One clear advantage of high-energy obser- 
vations is that photoelectric absorption has minimal ef- 
fects, so transmission-dominated CT AGN can be easily 
detected. It is only when the source becomes reflection- 
dominated that the emission at E>10 keV is affected. 

Current observations at E>10 keV are available only 
at relatively high fl uxes, and hence low redshifts, z<0.05. 
While BeppoSAX (jBoefla et all |1997[) was successfully 
used for targeted observations of known Seyfert galaxies, 
it is only now thanks to the International Gam ma-Ray 
Astro physics L aboratory (INTEGRAL; Winkle rTt al.l 
I2003f) and Swift (jGehrels et' al.l2004f) satellites that large- 
area surveys at these energies have been done. 

Usi ng the IBIS coded-mask telescope (jUbertini et al.l 
[200l . INTEGRAL surveyed ~80% of the sky down 
to a flux of 5 mCrab i n the 17-60 keV. The catalog 
of IKrivonos et al.l ()2007l ) reports the properties of 130 
sources detected in these all-sky observations and clas- 
sified as AGN. A large number of unidentified sources 
remain in this catalog, 48, but only seven are found a 
high galactic latitude (|6|>5°), and thus of likely extra- 
galactic origin. Five of the 130 AGN are CT AGN. 

We carried out a very deep survey with INTE- 
GRAL/IBIS, with a total exposure time of ~3 Msec, in 
the XMM-Large Scale Survey (XMM-LSS) region, reach- 
ing a flux limit of ~3xl0~^^ erg cm~^s~^ in the 20-40 
keV band (S. Virani in prep.). A total of 15 sources, all 
AGN, are found in this survey, including the prototyp- 
ical CT AGN, NGC 1068. We also found another CT 



AGN candidate, not detected in X-rays before. How- 
ever, an accurate Nh determination has not been ob- 
tained for this source yet, and hence it is not included 
in our sample. NGC 1 068 w as also included in the cata- 
log of I Krivonos et al] ()2007f ) and hence is already in our 
sample. 

Recently, iTueller et al.l (|2008( ) presented a catalog 
of 103 AGN detected in an all-sky survey with the 
Swift/BAT telescope. The 14 sources classified as blazars 
and BL Lac were excluded from our sample. 89 of the 
remaining AGN are at high galactic latitudes, |5|>15°, 
where only one source remains unidentified. The frac- 
tion of unidentified sources is much smaller for Swift 
compared to INTEGRAL because of follow-up obser- 
vations with Swift's dedicated X-ray telescope. In the 
ITueller et all (|200M ) catalog there are five AGN with esti- 
mated Nh greater than 10^^ cm~^. However, we caution 
that these Nh measurements were obtained by fitting a 
single absorbed power-law to the X-ray spectrum, while 
it is known that he avily absorbed AGN have more com- 
plex spectra (e.g., IVignati et al.l [19991 iLevenson et al.l 
[200l . 

We added to our sai nple the source NG C 7582, which 
has iV//'-1023 cm-2 in ITueller et"all (|2008D but has been 
shown with XMM/Newton data to have a very complex 
spectrum with strong evidence f or mildly Compton-thic k 
absorption, iV/z^lO^'^ cm'^ (jPiconcelli et all l2007f l. 
With the improved sensitivity of the Suzaku telescope 
(jMitsuda et al.l[2007f ). it is possible to perform detailed 
X-ray spectroscopy for some of the sources detected 
by Sw ift/BAT included in the catalog of Tuellcr et aT| 
(|2008f ). revealing in some cases Compton-thick absorp- 
tion. We added to our sample the sourc e NGC 5728, 
which as reported bv lComastri et al.1 (|2007t l from Sukazu 
observations is obscured by a Compton-thick gas with 
iV/f~2.1xl02'* cm"2. We also added the source ESO 
005-G004 wh ich base d on t he Suzaku observations re- 
ported by lUeda et al.l (|2007D is a heavily obscured, CT 
AGN. 

In summary, INTEGRAL and Swift found 130 and 
103 AGN, respectively, in their wide area surveys; 76 
sources (^58%) were detected by both surveys. (This 
fraction is not larger due to the differences in sky cover- 
age and the non-uniform depth of the observations.) We 
then found 15 AGN in deep 3 Msec INTEGRAL observa- 
tions, one of them the CT AGN NGC1068. INTEGRAL 
detected five CT AGN, while Swift found eight. How- 
ever, there is incomplete overlap between the two sam- 
ples and we note that the disparate energy ranges make 
direct comparison of fiuxes difficult. The INTEGRAL- 
detected CT sources are: NGC4945, Circinus, Markar- 
ian 3, NGC3281 and NGC 1068; the Swift/BAT CT 
AGN are: NGC4945, Markarian 3, NGC3281, NGC7582, 
NGC5728, NGC5252, NGC6240 and ESO 005-G004. 

The number of CT AGN found by these surveys is 
surprisingly low, compared to the sample of known CT 
AGN in the local Universe. In a study of optically- 
selected Seyfert 2 ga laxies with hard X-ray information, 
iRisahti et all (|1999l ) found 16 CT AGN in a total of 45 
Seyfert galaxies, although four were later demonstrated 
to most likely not be Compton t hick (NGC 1386, IC 
3639, NGC 5005 and NGC 4939: Maioli no et all [l998l : 
Ghosh et all 120071 : iGallo et al. ,2006) . Of the remaining 
12 CT sources, three were detected by Chandra and/or 



The Space Density of Compton Thick AGN 



3 



XMM, while the rest are mostly reflection-dominated 
sources, too faint to be detected by either INTEGRAL 
or Swift even tho ugh they are nearby , moderate luminos- 
ity AGN. In fact. lAwaki et all ()2008[ ) recently confirmed 
the CT nature of N GC 2273, one of the sources in the 
iRisaliti et al.l ()1999l ) sample, which is however too faint 
to be detec ted by INTEGRAL or Sw ift. 

Recently, iDella Ceca et al.l ()2008aD published a list of 
18 CT AGN with detections at E>10 keV. Of these, 
seven were detected by Swift and/or INTEGRAL, while 
the remaining 11 were studied with pointed BeppoSAX 
observations, and are typically fainter than the INTE- 
GRAL/ Swift detection threshold. We use all these sam- 
ples, suitably amended as necessary, to place constraints 
on the number density of CT AGN. 

An alternative way to find CT sources is by studying 
the water maser emission in AGN. Because large amounts 
of molecular gas are required to produce the maser am- 
plification, AGN with detected water maser emission are 
in g eneral heavily o bscured along the line of sight. In 
fact. lGreenhill et al . (2008) recently reported that from a 
sample of 42 AGN known to show water maser emission, 
95% have Nh>10'^^ cm'^ and 60% are Compton-thick. 
Since these AGN were not detected at high energies by 
either INTEGRAL or Swift we do not include them here; 
however, we note that water maser emission appears to 
be a highly efficient way to identify a large number of 
heavily obscured sources. 

2.1. The Log N-log S Distribution 

Figure [U shows the cumulative number counts of AGN, 
with CT sources shown separately, as a function of hard 
X-ray flux. In order to avoid the necessity of specifying 
a standard spectrum to convert fluxes to different energy 
bands, we show the INTEGRAL and Swift sources sep- 
arately, but note that a good agreement (within ^40%) 
in the normalization between the two distributions ex- 
ist if a standard band conversion is assumed. At these 
high fluxes the slope of the log N-log S is Euclidean, im- 
plying an uniform spatial distribution, as expected given 
the low rcdshifts of these sources. We also compare with 
the distribution predicte d by the AGN p o p ulation syn- 
thesis model with which iTreister fc Urrvl ()2005i) flt the 
XRB, and find in general good agreement in slope and 
normalization. 

The log N-log S relations for the five CT AGN de- 
tected by INTEGRAL and the eight sources observed 
by Swift/BAT are also consistent with Eculidean slopes, 
with normalizations of 10""' deg~^ at fluxes of ^5x10"^^^ 
and -9x10"" erg cm^^g-i in the INTEGRAL and 
Swift bands, respectively. (For a power law spec- 
trum with F = 1.8, for example, i^i4~i95,5'uii/t 
2.3Fn^Q0jNTEGRAL for most column densities Nh < 
24 cm"^.) About twice as many C T AGN at low r e dshift s 
were reported in the sample of IRisaliti et al.l (|1999( ). 
Since all these CT AGN have detections at high energies, 
we plot them in Fig. [T] Clearly, they fall significantly 
below the extrapolation of the observed log N-log S, sug- 
gestin g high levels of incompleteness in the IRisaliti et al.l 
()1999| ) sample. This is not surprising given that these 
sources do not come from a fiux-limited survey but from 
pointed observations. 

A possible source of incompleteness in our sample 
of CT AGN comes from the difficulty in measuring 
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Fig. 1. — LogN-logS distribution for AGN detected at high en- 
ergies. The gray line in the top panel shows the AG N in the well- 
defined Swift/BAT samples in the 14-195 keV band (ITueller et al.l 
[20M. while the bot tom panel shows the INTEGRAL sources 
IIKrivonos et al.||2007t ) in the 17-60 keV band. Solid squares show 
the 15 sources detected in the deep 3 Msec INTEGRAL observa- 
tions of the XMM-LSS field (S. Virani in prep.). Solid circles mark 
the CT AGN detected with Swift (top panel) and INTEGRAL 
(bottom panel). The black solid lines show the expected AGN 
logN-logS from the population synthesis model of Trcistcr & Urr^ 
poo s), which at these fluxes corresponds to a Euclidean distribu- 
tion. The dashed lines mark the Euclidean slope normalized to 
the number of Swift and INTEGRAL CT AGN. The gray lower 
limits show the previously-known transmission-dominated AGN 
with hard X-ray observations, not detected in the INTEGRAL or 
Swift surveys. These are lower limits since they were selected from 
pointed observations and are thus highly incomplete. 



the amount of absorption in these sources, given that 
they hav e in general very c omplex X-ray spectra. As 
noted bv ITueller et~al] (|2008D . sources that are not well- 
characterized in X-rays by an absorbed power-law are 
good candidates to be heavily obscured AGN. In their 
Swift/BAT sample, a total of 46 sources with complex 
spectra were reported. Of those, 18 have an optical clas- 
sification of Seyfert 1.5 or lower, and hence it is very 
unlikely that they are CT AGN. Considering the very 
extreme assumption that the remaining 28 sources are 
aU CT AGN increases the normalization of the CT AGN 
log N-log S by only a factor of ~2. This is because 
a large fraction of the complex-spectrum sources have 
fluxes fainter than that of NGC6240, one of the faintest 
confirmed CT AGN in the Swift sample. In any case, 
it is imp ortant to rern a rk th at according to a detailed 
study bv IWinter et al.l ()2009f ). while all these complex- 
spectrum sources are highly obscured, only half a dozen 
have some evidence of Compton-thick column densities. 
Hence, we conclude that the observed log N-log S is not 
affected significantly by this possible source of incom- 
pleteness. 

For the transmission-dominated AGN in our sam- 
ple (i.e., excluding NGC1068), we find volume densi- 
ties for Lx>1042 erg g-i 5.5t^:?xl0-5 Mpc'^; and 

for Lx>10'^^ erg s'^ of 2.2t?;?xl0-^ Mpc'^. Because 
this is a flux-limited sample, luminosity and redshift are 
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strongly correlated. For example, a source with X-ray 
luminosity of 10''^, 10''^ or 10'''* erg can only be de- 
tected up to z~0.005, 0.015 or 0.045, respectively. Thus, 
the source densities inferred here are valid only up to 
these limiting redshifts, corresponding to distances of 
~21, 63 and 190 Mpc. Also, because there is a sig- 
nificant correlation bet ween luminosity a n d fraction of 
obscu red sources (e.g., lUeda et al] 120031 : iBarger et al.l 
|2005[) . CT AGN are found only up to 0=0.024 in this 
sample, even though unobscured sources have been found 
by Swift/BAT up to z~0.1. That is, CT AGN preferen- 
tially have low luminosities, so they are found mostly at 
low redshift. In effect, the flux limit prevents us from 
detecting the many CT AGN at higher redshift. The 
derived volume densities for CT AGN at z:^0 are fully 
consiste nt with the values derived by iDella Ceca et al.l 
(|2008bD from three INTEGRAL sources only. 

Taking into account the densities reported here, 
the number of CT AGN relative to the X-ray- 
selected AGN population is 5.3x10-^2.2x10-^=24% 
for sources in the Lx=10'*^-'*'^ erg s-* range, while 
for sources with Lx = W^^~^'^ erg s-* this fraction is 
2.2x10-^2.9x10-^ =7.5%. This calculation uses the lu- 
minosity function of lUeda eF al. (200 ^. However, s imilar 
numbers are obtained if the lLa Franca et al.l (|2005l ) lumi- 
nosity function is used instead. Hence, the relative frac- 
tion of CT AGN decreases by about a factor of 3 for an 
order of magnitude increase i n luminosity. F o r Com pton- 
thin sources, according to iTreister et al.l ()2008[ ). the 
fraction of obscured sources {Nh > 10^^ cm-^) de- 
creases from 100% at Lx=10^^ erg s^^ to ^^60% at 
Lx—W"^^ erg s-*, implying a decrease of about a factor 
of 2. Therefore, the decrease in the fraction of CT AGN 
is comparable to the decrease in the fraction of obscured 
Compton-thin AGN, a reasonable agreement given the 
statistical errors in ou r sample. This is in agreement 
with the conclusions of iFiore et al.l ()2009f ). who found a 
similar decrease in CT AGN with increasing luminosity 
in their high-redshift IR-selected sample. 

2.2. Nh Distribution 

A key ingredient in our understanding of the AGN 
population and of the properties of the obscuring ma- 
terial is the distribution of line-of-sight column densities, 
parametrized in terms of the neutral hydrogen column 
density, Nh- In Fig. [5] we show the observed Nh dis- 
tribution for the s ources in the Swift/BAT sample of 
iTueller et aH (|2008f ) obtained from very simple spectral 
fitting assuming an intrinsic absorbed power-law spec- 
trum; the distribu tion from the AGN population synthe- 
sis models of , Trei ster &: Urrvl (j2005| ). adapted to the flux 
li mit of the Sw ift/BAT sample; the distribution assumed 
bv iGilli et al.l (|2007D : and the distribution pr e dicted by 
the galaxy evolution models of I Hopkins et al.l (|2006D . 

The Nh distribution observed in the Swift sample is 
relatively flat, before a strong decline at A^//>10^'' cm-^. 
This decline corresponds to a relatively low number 
of CT AGN in this sample, as mentioned before. In 
contrast, the Nh distribut ion used in the XRB pop- 
ulation synthesis model of ITreister fc Urrvl (|2005[ ) had 
roughly the same number of sources with Nh in the 
lO^'^-lO^^ cm^^ and 10^''-10^^ cm-^ ranges, because of 
an incorrect assumption about the normalization of the 
HEA01-A2 X-ray background. This translates into a 
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FlG. 2. — Distribution of neutral hydrogen column density (Nh) 
for the AGN detected in the Swift/BAT survey (solid histogram). 
The distribution is to first approximation flat, but shows a signif- 
icant decrease in the number of AGN with Nj{>lO^'^ cm~^. The 
dotted histogram shows the Nh distribution used in the IGilli et al.l 
II2007I) AGN population synthesis model, normalized to the num- 
ber of sources in the Swift/BAT survey, while the dot-dashed line 
shows the Nff distribution predicted by the galaxy- merger mod- 
els of IHopkins et al.l II2006I ). The dashed histo gram shows the 
N[{ distribution assumed by the ITreister fc Urrvl ( 2005) model for 
1020<Ar^<1024 cm-2, modified for the flux limit of the Swift/BAT 
survey and normalized to the same number of sources. The dis- 
crepancies at low Nh are not relevant to the present work (see 
§2.2 for details). For Nh>W'^'^ cm"^^ the dashed region shows the 
number of CT AGN allowed by the current observations. 



discrepancy of a factor of ~3 more CT AGN assumed by 
that model than is observed in the BAT sample. A simi- 
lar unrealistica lly high number of CT AGN was assumed 
in the work of IGilli et al! (|2007i ). The relationship be- 
tween the number of CT AGN and the XRB is explored 
in detail in §3 below. 

For relati vely unobscured sources, NH<W^^cm~'^, the 
iTueller eT^l . (2008) Nh distribution is significantly dif- 
ferent from that observed in the deepest X-ray fields. 
For example, for the sources in the Ch andra Deep Field 
North and South, ITreister et al.l (|2004f ) reported a sharp 
peak at 7V//=10^° cm-^ and almost no AGN in the 
10^°-10^* cm-^ range. The discrepancy at low values of 
Nh between the Swift/BAT and the deep fields samples 
could be due to the difficulty in measuring low Nh val- 
ues at higher redshifts, as discussed by e.g.. lAkvlas et all 
(|2006f ). The Nh distribution in the population synthe- 
sis model of [Treister fc Urrvl ([2005) matches well the 
observed distribution in the Chandra deep fields, but 
has some sign ificant differences with the distribution of 
ITueller et al.l ([2008). This discrepancy however is not 
relevant for our present work, which focuses on the most 
obscured AGN. 

2.3. Comparison with Models 

The Nh distribution fo r AGN is predicted by the 
galaxy evolution models of lHopkins et al.l (|2006[ ). They 
assumed that AGN are fueled solely by mergers of gas- 
rich galaxies, and the Nh distribution was derived by 
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integrating the amount of gas along the line of sight for 
the simulated galaxies. The resulting distribution peaks 
at Nh^IO^^'^ cm~^ and declines strongly towards higher 
column densities, in reasonable agreement with the num- 
ber of CT sources reported here. It is interesting to note 
that in the Hopkins et al. model the obscuring gas is 
located ~100-200 pc from the nucleus. Such a large 
scale for the obscuration disagrees with observations of 
a few nearby AGN using near-IR interferometry, which 
show outer radii for th e obscuring material of ~3 pc for 
NGG10 68 (iJaffe et al.l[20 04) and ~2 pc for the Circinus 
galaxy (|Tristram et al.[200'it ). for example. Similarly, the 
latest torus models predict small scales fo r the obscuring 
material, <10 pc (|Nenkova et all l2008bD . although fit- 
ting the IR AGN spectra does not provide a very strong 
constraint to the torus size. Further comparisons will 
provide an interesting test of the Hopkins et al. model. 

The fact that a relatively small number of AGN with 
A^//>10^'* cm~^ is observed can be used to constrain the 
nature of the obscuring material. This lack of CT AGN 
can b e interpreted either i n the conte xt of a clumpy to rus 
(e.g. jKrolik fc Begelma^ [l988: Nen kova et al.ll2008al) or 
a smooth distribution (e.g., iPier fc^rolikl 19921). Fo r 
example, the Nh distribution of iTreister fc Urrvl (|2005D . 
presented in Fig.|2l assumed a smooth torus with a single 
equatorial column density of lO^'^ cm~^. To accommo- 
date a smaller number of CT AGN while still matching 
the observed Nh distribution for Compton-thin sources 
would require a distribution of equatorial densities, in 
which only a small fraction of the AGN reach the CT 
levels for nearly-equatorial line of sights. In the case of a 
clumpy torus, the explanation is perhaps more natural; 
the small fraction of CT AGN implies that only a few 
sources have a la rge number of clouds, e .g., > 10 clouds 
for the models of INenkova et al.l (|2008af) . 

3. CT AGN AND THE X-RAY BACKGROUND 

3.1. Parameter Degeneracies 

The spectrum of CT AGN at high energies, even 
for transmission-dominated sources, is often domi- 
nated by the Compton reflection component (e.g., 
IMatt et al.ll2000l). which has a strong peak at E~30 keV 
(jMagdziarz fc Zdziarski|[r995D . The observed spectrum 
of the XRB, which we now know is just the integrated 
emission from previously unresolv ed AGN, also has a 
peak at about the same energy (|Gruber et al.l I1999D . 
Hence, it was suspected for a long time that CT AGN 
provide a significant contrib ution to the XRB emis sion. 
In fact, in the early work of IComastri et all (|1995[ ) the 
contribution of CT AGN to the XRB is -20%, similar 
to the value assumed in t he popu lation synthesis mod- 
els of fUeda et all fl2003h. ITreister fc Urrvl (|200,^ . and 
iGiin et all (|2007f) : IShankar et al.l (|2009D report a shghtly 
higher contribution of ^^30% at ~20 keV. 

Because it is very hard to measure the number den- 
sity of GT AGN, even locally, AGN population synthesis 
models have assumed it to be a fixed fraction of the ob- 
scured, Compton-thin sources, typically ~0.5-l times as 
many. In Figure |3| we show the fraction of all AGN that 
are Gompton-thick, compared to the observed fraction in 
the INTEGRAL and Swift samples as a function of hard 
X-ray flux. At fluxes of ^10 ~^^ erg cm~^ ~^ th e fraction 
of CT AGN in the model of iGilh et al] (pOOl is ~15%, 
while the observed value is ~6±5% for INTEGRAL, and 
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Fig. 3. — Measured fraction of AGN that are Compton-thick in 
the INTEGR AL ( open circle , Krivonos ct al. 2007) and Swift/BAT 
{filled circles; ITueller et al.l 12008.) samples. The filled circle with 
solid error bars shows the fraction using the identified sources only, 
while the circle with dotted error bars assumes that the one uniden- 
tified source is a CT AGN. The solid line sho ws the fraction of 
CT AGN from the modified ITreister fc Urrvl pOOS. ) population 
synthesis model; the original assumption was a factor of ~4 too 
high, so was modified as described in the present text to match the 
Swift/BAT and INTEGRAL observation s. The da shed line shows 
the fraction of CT AGN in the model of lGiUi et a l. (2007), which 
is a factor of ~3 higher than observations, and increases sharply 
at faint fluxes because of the assumed steep dependence of the 
fraction of obscured sources on luminosity. 



~8±3% for the Swift sample. For comparison, in Fig- 
ure [3] (solid line) we show the predicte d CT AGN frac- 
tion a s a function of flux for the model of lTreister fc Urrvl 
()2005f l with the number of GT AGN modified to match 
the INTEGRAL and Swift observations presented here. 
The lTreister fc Urrvl (|2005l ) model assumes a nearly con- 
stant fraction of GT AGN, while the iGilh et all |2007l ) 
model assumes a steep increase in the number of CT 
AGN at fluxes fainter than ~10~"'^'' erg cm~^ s'-^. This is 
due to the assumed lumin osity dependence of the fraction 
of obscured sources in the lGilli et al.l (|2007l ) model, which 
decreases steeply above luminosities of ~10^^ erg s^^ and 
is flat at lower luminosities. While such faint fluxes are 
still out of reach for current hard X-ray observatories, it 
will be possible to test this flux regime with NuSTAR and 
the International X-ray Observatory (IXO). However, the 
luminosity dependence of the fraction of obscured AGN 
assumed by the model o fiGilli et al.l (|2007f ) can already be 
ruled out, in particular at high luminosities, b y observa- 
tions of Compton-thin AG N at lower energies (jHasingeil 
l20f)8l : lTVeister et al.ll2008l) . 

The fact that the XRB does not constrain the num- 
ber density of GT AGN can be explained by strong de- 
generacies in othe r parameters, like the assumed sprea. d 
in spectral index (jGiUi et al.ll2007HShankar et al.ll2Ci09h . 
high-energy cutoff, etc. When trying to constrain the 
number density of CT AGN the most important pa- 
rameter is the normalization of the Compton reflection 
component, which is directly related to the AGN lu- 
minosity at ^^30 keV, where the CT AGN contribu- 
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tion is maximal. Given that even the nearest AGN 
have only low signal-to-noise observations at E>10 keV, 
the normalization of the Compton reflection component 
is not well constrained by observations of individual 
sources. Fro m a sample of 22 Seyfcrt galaxies, excluding 
CT sources, iMalizia et al.l ()2003.) concluded that both 
obscured and unobscured sources have similar reflec- 
tion components with normalization values in the range 
R'^O.G-l (in units of 2tt). A sim ilar value of R~l was 
reported by iPerola et al.l (|2002D based on BeppoSAX 
observations of a sample of 9 Seyfert 1 galaxies. Al- 
though with large scatter, normalizations for the aver- 
age reflection component of 0.9 for Seyfert 1 and 1.5 
for Seyfert 2 were measured from BeppoSAX observa- 
tions of a sample of 36 sources (jPeluit fc Courvoisieil 
l2003f ). Early population synthesis models for the XRB 
assumed values of R=1 .29 for unobscured sources and 
0.88 for obscured AGN (jComastri et al.lll995l: i GiUi et al 
19991). In contrast, t he m odels of lUeda eFa l. (200 ^^ 



Treister fc ([2003) and iBallantvne et al.l (|2006) as 



sumed a constant value of R=l (equivalent to a solid 
angle of 27 r ) for b oth obscured and unobscured sources. 
Gilli et aTl il2007tl assu med the same normalizations as 
Comastri et al.l ()1995D . R'^1.3 and 0.88; however, for 
high-luminosity sources, ix>10'*^ erg s^^ , the reflec- 
tion component was not included (R=0). 

For a given number of CT AGN, the resulting inten- 
sity of the XRB at ~30 keV is directly linked to the 
assumed normalization of the reflection component. In 
Figure |4] we show the values of the normalization of 
the Compton reflection component and CT AGN num- 
ber density that produce XRB intensities in the 20- 
40 keV region consiste nt with the latest ob served val- 
ues from INTEGRAL (IChurazov et al.l[200l and Swift 
(|Aiello et al.ll2008f) . This is the energy range in which 
the contribution of CT AGN to the XRB is maximal 
and hence can be best constrained. F or comparis o n, the 
parameters assumed by the model of iGilli et al.l ()2007f ) 
lie in the upper left region of Figure [H at a density of 
CT AGN roughly 3x higher than the observed value of 
~2xl0^^ Mpc~^ and average Compton reflection com- 
ponent normalization of ~0.6. (The latter is a rough 
estimate, since they assumed R = 0.88 for obscured 
sources at low-luminos ities, and R = at h igh lumi- 
nosities.) The model of iTreister fc Urrvl (|2005 f) assumed 
a similarly high number of CT AGN and a higher nor- 
malization of the Compton reflection component, and 
hence resulted in a highe r intensity of the X RB, con- 
sistent with the HEAO-1 (jCruber et al.lll999D measure- 
ments increased b y 40%, si milar to what was assumed by 
lUeda et al.l (|2003D and Ball antvne et all (2006) in order 
to match the observations of the XRB at lower energies 
by Chandra and XMM. Such a high value of the XRB 
intensit y at E^lO-50 keV is no w ruled out by new INTE- 
GRAL (IChurazov et al.ll200l. Swift JAiello et al.ll20Ml 
and BeppoSAX (jFrontera et all I2007D data. Given the 
degeneracies with other model parameters, it is unlikely 
that the XRB could be used to constrain the average 
value of R. High signal-to-noise observations of individ- 
ual sources at E>10 keV are required for this purpose. 

3.2. A New X-ray Background Fit 

Since both the number density of CT AGN and the 
normalization of the Compton reflection component can 
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Fig. 4. — Degeneracy of the local density of CT AGN and the 
normalization of the Compton reflection component subject to the 
constraint of either the XRB intensity or the number of CT AGN in 
hard X-ray surveys. The dark gray region shows the space density 
obtained from the 10 CT AGN detected by Swift and INTEGRAL 
in complete samples at z~0, including l-a statistical fluctuations. 
The light gray region shows the constraints to these parameters 
given by the intensity of the XRB in the 20-40 keV band, assuming 
a 5% uncertainty in the XRB intensity. The correct values of these 
parameters must be at the intersection between these two regions, 
namely, a normalization of the Compton reflection component of 
~1 and a CT number density of ~2xl0~^ Mpc""^. 



now be constrained independently, we can attempt to 
match the observed spectrum and intensity of the XRB. 
In Figure [5l we show our new fit, which matches the 
INTEGRAL and Swift observations at E>10 keV, which 
are ~10% hi gher than the HEAO-1 normalization. The 
original fit of ITreister fc Urrvl (|2005f) . which has a factor 
of ~4 more CT AGN, is also shown. Not surprisingly, 
the effects of changing the number of CT AGN are most 
important in the E=10-100 keV region. 

The new XRB fit matches both the INTEGRAL and 
Swift/BAT observations at E>10 keV and the Chandra 
measurements at lower energies (which are ^30% higher 
than the HEAO-1 A2 observations). Recently, a new 
measurement of the XRB intensity at E=1.5-7 keV us- 
ing the Swift XRT (X-ray telescope) was presented by 
iMoretti et aT] (|2009D . These new data confirmed that 
the original HEAO-1 normalization should be increased 
by ~30% and ~10% at low and high energies respec- 
ti vely. In con t rast, t he AGN population synthesis model 
of iGilh et"aI1 (|2007f ) assumed the original HEAO-1 in- 
tensity at all energies, which translates into a relatively 
lower contribution from unobscured so urces. In o r der to 
produce the necessary hard spectrum, I Gilli et al] ()2007f ) 
had to assume a relatively high number of obscured 
sources at high luminosities, i.e., an unusual, inverted 
dependence of the obscured fraction of AGN as a fun c- 
tion of luminosity ()Hasingeij|2008t ITreister et al]|2008l ). 

Assuming a fixed value of the Compton reflection com- 
ponent, how much can the number of CT AGN change 
and still match the XRB, given the existing uncertain- 
ties in the intensity measurements? The INTEGRAL 
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Fig. 5. — Observed spect rum of the extragalactic X- 
ray backsround from HEA O-1 dGruber et al.l l lOgg). C handr a 
lIHickox fc M arkevitch' '200^, XMM l IDe Luca fc Molcndi "200?), 
INTEGRAL (Churazov ct al. 2007) and Swift jAicllo ct al. 2008) 
data. The dashed gray line shows the XRB s pectrum from the 
AGN population synthesis model of lTreister &: U rrv (2005), which 
assumed a 40% higher value for the HEAO-1 XRB normaliza- 
tion. The thick black solid line shows our new population synthesis 
model for the XRB spectrum; the only change is the number of CT 
AGN, which is reduced by a factor of 4 relative to the number in 
ITre istcr & Urry (2005). Red, blue and thin black solid lines show 
the contribution to this model from unobscured, obscured Comp- 
ton thin and CT AGN respectively. 



measurements of the XRB, reported by IChurazov et ahl 
([2007), have uncertainties of ~5% including both statis- 
tic and systematic effects. Similarly, the Swift measure - 
ments have estimated errors of ~3% (lAiello et al.ll2008l) . 
Both measurements are consistent with each other but 
are ~10% higher than the original HEAO-1 intensity. 
Translating this '^5% uncertainty into an uncertainty in 
the number of CT AGN, we conclude that the total num- 
ber of CT AGN can be changed by a factor of 55% and 
still be consistent with the current measurements of the 
XRB. However, this calculation does not include the un- 
certainty in the normalization of the Compton reflection 
component, which is by far the dominant factor. For 
comparison, the statistical errors for the measurement of 
the ten CT AGN detected combining the Swift and IN- 
TEGRAL surveys correspond to an uncertainty of ^--^30% 
(jGehre ls 1986), i.e., the direct detection of CT AGN is 
much better than the XRB for determining the number 
of CT AGN. 

Given that the number of CT AGN in the local Uni- 
verse is effectively constrained by the Swift and INTE- 
GRAL surveys, it is now possible to estimate the total 
contribution of CT AGN to the XRB, as well as its red- 
shift dependence. In order to do that, we integrate the 
total X-ray emission from the CT AGN in our popula- 
tion synthesis model and divide it by the observed XRB 
intensity. To facilitate the comparison with the local 
sources observed by Swift, and to make sure that the 
effects of absorption are negligible, we perform this inte- 
gration over the 14-195 keV band. In Figure [6] we show 




Fig. 6. — Cumulative fractional contribution of CT AGN 
to the XRB in the 14-195 keV Swift/BAT band as a func- 
tion of redshift, from the population synthesis model pre- 
sented here {solid line; see text for details). As shown by 
the vertical dashed lines, 50%, 80% and 90% of the total CT 
AGN contribution come from sources at z<Q.7, 1.4 and 2.0 
respectively. Only ~1% of the total XRB intensity comes 
from CT AGN at z>2. Given the current 5% uncertainties in 
the measurement of the XRB intensity, this means that the 
XRB spectrum does not constrain the number of high-redshift 
CT AGN at all (factor of 5 uncertainty). The data point at 
zr^O corresponds to the contribution to the XRB by the CT 
AGN detected by Swift/BAT, while the data points at high 
redshift were obtained from the CT AGN in the Chandra 
sample of iTozzi et all ||2006D . Solid error bars correspond to 
transmission-dominated sources only, while the data points 
with dashed error bars include all the sources in the sam- 
ple. As expected, most of the contribution to the XRB comes 
from the transmission-dominated sources, which are in gen- 
eral brighter. Good agreement between our population syn- 
thesis model and observations of CT sources is found at all 
redshifts. 



the resulting redshift dependence of the fractional contri- 
bution to the hard XRB radiation. As can be seen, the 
total contribution of CT AGN to the XRB is --9%, and 
about 50% of it comes from sources at 2;<0.7. Similarly, 
we conclude that ~2% of the XRB is provided by CT 
AGN at z>1.4, while CT AGN at z>2 only contribute 
<1% to the XRB. Conversely, the 5% uncertainty in the 
absolute measurement of the XRB intensity translates 
into an uncertainty of a factor of ~5 in the number of 
CT AGN at z>2. Hence, the number of CT AGN at 
high redshift is largely unconstrained by the XRB. 

In Figure [6] we further compare this expected redshift 
dependence to the integrated fluxes from individually- 
detected CT AGN. At 2:~0 we integrate the emis- 
sion from the eight sources detected by Swift/BAT. At 
higher redshifts we use the sample of CT AGN can- 
didates detected in X-rays in t h e Cha ndra Deep Field 
South reported by iTozzi et al.1 (|2006[ ). which includes 
14 reflection-dominated AGN and six transmission- 
domi nated AGN with Nh> 1 0^'^ cm ~^. In the same 
field, iGeorgantopoulos et al.l ()2007D found a total of 
18 CT AGN candidates, but only eight of them with 
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a measured Nh greater than 10^^ cm^^; the remain- 
ing sources were selected based on their flat X-ray 
spectra. All s ix transmiss i on-do minated CT AGN in 
the sam ple of iTozzi et alj (120061) are included in the 
work of iGeorgantopoulos et al.l ( 20071 ). However, no 
overlap is found between th e reflection-domin ated CT 
AGN candidates reporte d bvlTozzi et al.1 (I2006D a nd th e 
flat-spectrum sources of IGeorgantopoulos et al.l (|2007D . 
Hence, it is possible that either selection of heavily ob- 
scured sources is highly incomplete. In order to compare 
properly with the local sample, in Figure[6]we show sepa- 
rately the contribution from the transmissio n-dominated 
AGN and from all sources in the sample of iTozzi et al.l 
(|2006| ). We separated the sample at z=1.5, to have 
the same number of sources in each redshift bin. As 
expected, most of the contribution to the XRB comes 
from the transmission-dominated sources, which are in 
general brighter in the X-ray band. The agreement at 
low redshift is not surprising, since by construction our 
model was adjusted to match the Swift/BAT observa- 
tions. However, it is very interesting that also for the 
high redshift sources the calculated contribution of CT 
AGN to the XRB agrees well with the limits from deep 
surveys. 

4. HIGH-REDSHIFT CT AGN 

As shown in the previous section, the number of CT 
AGN at high redshift is largely unconstrained by the 
XRB or by current hard X-ray surveys. Since a large 
fraction of the absorbed energy in heavily obscured AGN 
is re-emitted at mid-IR wavelengths, deep Spitzer data 
have been used to flnd CT AGN candidates at z>2. 
iDaddi et al.] (|2007[ ) used the excess in mid-IR luminosity 
(compared to UV estimates of star formation rates) to 
find obsc ured AGN n ot individually detected in X-rays. 
Similarly, iFiore et al.l ((2008.) used a combination of red 
optical-to-near-IR colors and high 24 /im luminosity to 
select CT AGN candidates at z~2. In both cases, very 
high source densities ha ve been estimated for mid-IR- 
selected CT AGN, e.g., iDaddi et all (|2007D reported a 
sky density of ~3,200 deg~^, similar to that of all pre- 
viously known AGN at those redshifts in the Chandra 
deep fields. Somewhat surprisingly, very little overlap 
is found between the two selection methods, implying 
the possible presence of interlopers. If these candidates 
are confirmed, a very large number of CT AGN exist at 
high redshift, considerably larger than the local popula- 
tion. This is qualitatively c onsistent with th e evolu tion 
in obscuration detected by iTreister &: Urrvl (l2006l) for 
Compton-thin sources. Recently. I Alexander et al! nt2008h 
reported the confirmation of seven CT AGN using opti- 
cal and mid-IR spectroscopy in the Chandra Deep Field 
North region, implying a similarly high density for CT 
AGN at high redshift. 

In order to quantify the density of CT AGN at high 
redshift, and to compare with local observations and 
AGN luminosity functions, in Figure [7| we present the 
available measurements of the comoving volume den- 
sity of CT AGN candidates as a function of redshift. 
Our measurement of the density of CT AGN at z—Q is 
~2.2xl0~^ Mpc~^, for sources with Lx >10'^^ erg s^\ 
as shown in §2.1. At higher redshifts and luminosi- 
ties, we infer the space density from several samples. 
Five X-ray-selected CT AGN candidates were found 
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Fig. 7. — Comoving density of CT AGN as a function of redshift 
in several luminosity b ins. Measured values (details in text) are 
shown by: open circle: [PoUetta ct al.' (^006), squares: 'Tozzi et al 



(2006), filled triangle: [P addi ct al. 
12009). star: (Alexander ct al. 



^ 



(2007), pentagons: Fiorc ct al 
(2008), filled circle: Swift/BAT and 
INTEGRAL, for luminosity limits indicated by color: magenta: 
Lx>10**5 erg s-\ blue: Lx>10''*' erg s-\ red: L^>10''3 erg s"! 
and green: Lx>10'*^ erg . Solid line s show the space densi- 
ties at corresponding luminosities for the ITreister fc Urrvl ()2005f ) 
model including the evolution of obscured AGN reported by 
ITreister fc Urrv) l (2"006, ). 



with the number of CT AGN adjusted to 
the observed local value (present paper). At high redshifts, a rel- 
atively large density of CT AGN is observed, compared to expec- 
tations from the evolving luminosity function. This suggests that 
the local sample is incomplete; or that the high-redshift infrared- 
selected samples include a large number of interlopers; or that CT 
AGN follow a different evolution than Compton-thin sources. 



in the Spi tzer Wide-area Infra Red Extragalactic survey 
f SWIRE: IPolletta et al.l 12001 : these are transmission- 
dominated CT AGN and all of them have intrinsic X-ray 
luminosities greater than 10"*^ erg s~^ (magenta open cir- 
cle . Fig. III). Twenty X -ray-selected CT AGN were found 
bv lTozzi et all ()2006[) in the Chandra Deep Field South; 
in this case, we separated the sample into two redshift 
bins and computed the comoving number density sepa- 
rately for sources with Lx>10'*'^ and Lx>^^'^'^ erg s~^ 
(squares. Fig. [7]). Finally, we also show the comoving 
number density estimate d from the mid-IR CT AGN can- 
didate s in t he samples of IDaddi et al.l (|2007l ). lFiore et all 
(|2009f) and I Alexander et al.l (|2008[ ). 

The expected comoving number density for CT AGN 
as a function of redsh i ft wa s computed using the mod- 
els of ITreister fc Urrvl ()2005f ). Briefly, we used the hard 
(2-10 keV) X-ray l umino sity function and AGN evolu- 
tion of lUeda et al.l (|2003f ) , with the Njj distribution and 
lumi nosity dependence of the fraction of obscured AGN 
of ITreister fc Urrvl (|2005D . In addition, we include the 
evolution of the relative nu mber of obscured sources re- 
ported bv ITreister fc Urrvl (|2.0G6). The normalization of 
the relative number of CT AGN was chosen to match the 
observed numbers at Z'^O. This model fits the observed 
XRB spectral shape and normalization, as shown in Fig- 
ure |5| The resulting comoving density of CT AGN as a 
function of redshift is shown in Figure |7| for sources with 
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Lx>10'^^ lO'^^, and lO'^'* erg s^\ 

While the observed density of X-ray-selected CT AGN 
with Lx>10^'^ erg s~^ is in pretty good agreement with 
the expectations at all redshifts, at higher luminosities 
the observed values are mostly higher than the expecta- 
tions. In fact, the comoving density for Lx>W^^ erg s~^ 
sources from SWIRE is almost two orders of magnitude 
higher than the expected value. However, it is impor- 
tant to note that this comoving density is derived from 
CT AGN candidates based on the spectral properties de- 
rived from low signal-to-noise observations, and also the 
number of sources detected is small, so the uncertain- 
ties are large. The values for X-ray-selected CT AGN 
with Lx>W^'^ erg s~^ are also higher than expectations, 
but in this case only by factors of ■^2-3. Similarly, the 
densities inferred from mid-IR-selected CT AGN are sys- 
tematically higher than the expected values at all lumi- 
nosities, typically by one order of magnitude. 

This discrepancy between expectations and observa- 
tions at high redshift can be interpreted in several ways. 
One obvious possibility is that the observations at low 
redshift are missing a significant number of CT AGN, 
which are included in the high-redshift samples. In fact, 
we have shown before that the high-energy surveys per- 
formed by INTEGRAL and Swift are mostly complete 
for transmission-dominated sources, but miss a signifi- 
cant fract ion of the reflection-d ominated AGN. While the 
sample of iRisaliti et al.l (|1999[ ) includes these sources, it 
is based on pointed observations and hence it is highly 
incomplete as well (Fig. [1]). Also, it is very likely that 
the high-redshift IR-sclccted samples include both trans- 
mission and reflection-dominated sources. Another pos- 
sibility is that the high-redshift samples, both X-ray- and 
mid-IR-selected, include a significant number of interlop- 
ers. These could be either less obscured AGN in the case 
of X-ray selection, or non-active galaxies undergoing sig- 
nificant but dusty star formation, and thus showing high 
mid-IR luminosities not due to AGN activity. Finally, it 
is possible that CT AGN follow a different evolution than 
Compton-thin sources. It is important to note that the 
evolu tion in the ratio of obscu red to unobscured 
AGN found bv lTreister fc Urrvl (|2006[ ) for Compton-thin 
sources is already included in the predicted volume den- 
sities. Hence, if both the low and high redshift observed 
samples are not systematically missing a significant num- 
ber of sources, the excess of CT AGN at high redshift 
implies a different (stronger) evolution for these heavily 
obscured sources than for obscured but Compton-thin 
AGN. 

5. THE DENSITY OF SUPERMASSIVE BLACK HOLES 

Because AGN are powered by accretion of gas onto a 
supermassive black hole, the AGN lumin osity functio n 
represents the history of cosmic accretion (|Soltanlll98d) . 
Hence, the AGN bolometric luminosity can be converted 
into a mass accretion rate, assuming an efficiency for the 
conversion e=L/rhc^ (typically, eciO.l). Then, the co- 
moving black hole mass density can b e written (following 
equation 17 of lYu fc Tremainal2002t l as: 

(1) 



where ^(L,z) is the evolving AGN luminosity function, 
BC{Lx) is the bolometric correction starting from the 
2-10 keV luminosity, and f{NH,L) is the ''Nh func- 
tion," or the fraction of sources at a given luminosity 
with a given Nh- For this calculat ion we use d the 2- 
10 keV luminosity function of .Ueda et al.l (|2003l ) and the 
Nh function with a luminosity dependence described in 
Section 3.2 of iTreister fc Urrvl ()2005D . The bolometric 
correction was calculated using the spectral energy dis- 
tri bution of a complete ly unobscured AGN as specified 
bv ITreister et al.l (|2006( ). as appropriate for the unified 
model of AGN. Additionally, we updated the spectral 
library with the new relation between X-ray (at 2 keV) 
and UV (25GGA ) luminosities using the val ue of the slope 
of the power-law extrapolation reported bv lSteffen et al.l 

dlool), 

aox = (-0.077 ± 0.015) log(L2 keV) + (0-492 ± 0.387). 

With these assumptions, the bolometric correction 
ranges from ~25 for Lx—W^^ ergs s~^ to ^^100 for 
Lx=10^^ ergs s~-^, in approximate (factor of ~2) 
agreer nent with the values assumed by other authors 
(e.g., iKuraszkiewicz ct al.' '2003'; 'Marc oni et"al] 12004 
iBarger et al. 2005; Hopkins et al. 2007). 

In previous works (e.g., lYu fc Tremai^ l2002f) . the 
AGN luminosity function was integrated from Lmin=0 
to Lmax—oo, which leads to very large extrapolations 
in particular at the faint end. In the present case, we 
use th e same integration limits used bv lTreister fc Urrvl 
12005) in their AGN population synthesis model, 
namely, i„,„=10''^-^ ergs s"\ Lmax=W^^ ergs s"\ 
A^ff,rmn=10^° cm"^ and Nn.max^W^^ cm"^. The num- 
ber of CT AGN in this model is matched to the IN- 
TEGRAL and Swift results, as reported above. With 
these assumptions, and using the typical value of e=0.1, 
we obtain a value for the local black hole mass den- 
sity of p(z=0)=4.5x 10^ Mq Mpc-^ This calcula- 
tion agrees well with the values est imated from observa - 
tions: p=4.6ti;4xlO^ MqMpc"^ ()Marconi et"aI1 12004( ) 
and p=(3.2-5.4jxl05 MoMpc-^ (IShankar et al.|[2009l ). 

In Figure |8| we present the black hole mass density as 
a function of redshift estimated fro m our calculation, to - 
geth er with the curves presented bv lMarconT "eTall (f200i ) 
and [Ya_&_Trcmainc ( 2002) . The main differences with 
the work of IMarconi et al.l (|2004l ) are in the number of 
CT AGN (they assumed 4 times more), the assumed 
bolometric correction (our is ^^3 times higher at low lu- 
minosities) and the redshift limit of the integration. Of 
these, the bolometric correction dominates, such that our 
derived local black holes mass density is slightly larger, 
even with the reduced nu mber of CT AGN. Not e that the 
bolometric correction of IMarconi et al.l (|2004f ) was ob- 
tained from observations of high-luminosity sources only, 
while our bolometric correction was tested by obse rva- 
tions of fainter sources as well (iTreister et al.lll2006D . A 
remarkably good agree ment is found betwee n our results 
an d the recent wor k of IShankar et all ()2009l ). Compared 
to lYu fc Tremaind (|2002D . we find twice the local inte- 
grated black hole mass density, because they used an 
ojjtigjjiijr aj]g,sar luminosity function and evolution, which 
p'eaKS at a ' high er redshift, z ~2, and evolves strongly. 
In contrast, the lUeda et al.l (|2003f ) luminosity function. 
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Fig. 8. — Black hole mass density as a function of redshift (inset 
shows the same curves normalized to their values at z=0), assuming 
an efficiency e=L/mc^= 0.1. The solid lines show the evolution for 
the population synthesis model described in this p aper, while the 
dotted lines show a similar calculation presented bvlMarconi e t aLl 
m04). In both cases the AGN luminosity function of Ueda ct al] 
ll2003i) was used, and the only differences are the number of CT 
AGN (4 times more in the model of Marconi ct al. 2004) and the as- 
sumed bolometric correction (~ 3 times higher for our calculation). 
The gray rectangle at 2=0 sho ws the ran ge of values consistent with 
observations, as reported bv IShankar et al. ( 2009). For compari- 
so n, the dashed lin e s show the black hole mass density estimated 
by lYu &: Tremainj 1)20021 ') which considered only high-luminosity 
unobscured sources. 



which includes lower luminosity and obscured sources, 
peaks at z ~1.1. In our calculation, the vast majority 
of the black hole growth occurs at low redshift (~5G% 
from z = 1.3 to 0), whi ch matches observa tions of AGN 
detected in X-rays (e.g, iBarger et al1l200lD . 

The space density of CT AGN is consistent with, but 
cannot be constrained by, the observed local black hole 
mass density. In addition, analogous to the weakness of 
the XRB integral constraint on the number of CT AGN, 
numerous degeneracies with other parameters, like the 
assumed bolometric correction and efficiency, are impor- 
tant. Even taking into account only the uncertainties in 
the local black hole mass density, we could still increase 
the number of CT AGN in the local Universe by factors 
of ~3. Hence, we can conclude that direct observations 
of CT AGN at high energies, like the INTEGRAL and 
Swift observations discussed here, are currently the only 
way to constrain the population of heavily obscured su- 
permassive black holes. 

6. CONCLUSIONS 

In this paper we constrain the space density of CT 
AGN in the local Universe using the recently-available 
wide-area surveys at high energies performed by INTE- 
GRAL/IBIS and Swift /BAT. A total of ten CT AGN 
at z<0.03 were found by either INTEGRAL and/or 
Swift. These observations are complete for transmission- 
dominated CT AGN, but arc probably still missing heav- 
ily obscured sources with Nh>10^^ cm^^. We find that 
the space density of local CT AGN follow a Euclidean 



distribution with a normalization of ^^10 * deg ^ at 
fluxes of -5x10"" and -9x10"" erg cm ^s ^ in the 
17-60 keV and 14-195 keV bands, respectively. This is 
about 3-4 times smaller than the values expected from re- 
cent AGN population sy nthesis model s that flt the extra- 
galactic XRB ( Treister fc Urrvl [20051 : iGiUi et all [20071 
and thus modifications to those models is required. 

We present here a new population synthesis model for 
the XRB, with the number of CT AGN constrained by 
the INTEGRAL and Swift number counts. We find 
that the fraction of AGN that arc Compton thick at 
-p20-ioofcey~10~^^ erg cm~^s"^ is —5%. We show that 
the XRB by itself cannot be used to constrain the num- 
ber of CT AGN, mainly due to degeneracies with other 
parameters, the most important of which is the normal- 
ization of the Compton reflection component. We find 
that the total contribution of CT AGN to the XRB is 
-9%, with only <1% from CT AGN at z>2. Hence, 
taking into account the 5% uncertainty in the XRB in- 
tensity measurements, the number of CT AGN at high 
redshift is essentially unconstrained by the XRB, even if 
all the other parameters could be fixed. 

We calculate the local black hole mass density i nferred 
from AGN activity using Soltan's argument (jSoltanl 
|1982| ). taking into account the contribution from CT 
AGN estimated in this work. For an accretion efficiency 
e=L/rnc^= 0.1, we find an integrated local black hole 
mass density of 4.5 x 10^ Mq Mpc""^, in excellent agree- 
ment with recent estimates based on measured masses of 
local dormant black holes. Considering the current un- 
certainties in these estimates, we conclude that only the 
direct observations of CT AGN such as those discussed in 
this paper can effectively constrain the number of heav- 
ily obscured AGN. Based on the number density of CT 
AGN presented here, our best estimate of the fractional 
contribution of CT AGN to the total accreted black hole 
mass is < 10%. 

Using a combination of X-ray and mid-IR selection, 
the space density of CT AGN at high redshift is starting 
to be constrained. We find that observed densities are 
systematically higher than expected from the evolving 
AGN luminosity function measured from less obscured 
sources, assuming iV/f -independent evolution of the lo- 
cal CT AGN population. This can be explained in three 
ways, any or all of which could be the case. First, the 
local sample might be incomplete, particularly because 
even hard X-ray selection is biased against reflection- 
dominated CT AGN. Second, the high-redshift samples 
may be contaminated by strongly star-forming galaxies 
or other interlopers. Third, CT AGN may evolve more 
strongly than less-absorbed sources, implying a relatively 
larger number of CT AGN in the early Universe. To de- 
cide this question requires the help of observations with 
the new generation space-based hard X-ray observato- 
ries. While mid-IR selection of heavily obscured AGN 
is very promising, these samples are inevitably affected 
by the presence of interlopers, in particular from star- 
forming galaxies, and by the lack of accurate measure- 
ments of the amount of obscuration. 

Hard X-ray selection provides a cleaner sample of 
CT AGN, since Nh values can be measured directly 
and there is almost no contamination from star-forming 
galaxies at these energies. Several different approaches 
are currently being planned for the next generation of 
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high-energy (E>10 keV) missions, to provide a large and 
complete sample of CT AGN up to z^3. The Energetic 
X-ray Imaging Survey Telescope^ (EXIST; iGrindlavl 
l2005f ) will perform an all-sky survey in the 20-80 keV 
energy band to flux limits of ~6xl0~^^ erg cm~^s~^, 
finding thousands of heavily obscured AGN up to z~l 
and high-luminosity CT quasars at all redshifts. With 
a complementary approach, t he Nuclear Spe ctroscopic 
Telescope Array*^ (NuSTAR; iHarrisonl I2008D . with a 
scheduled launch date of August 2011, will perform tar- 
geted observations of fields of ~1 deg^ to flux limits of 
^2x10"^"* erg cm~^s~^, hence ^20 times deeper than 
EXIST, in the 6-79 keV band, for exposure times of 
~1 Msec; these observations will be able to find low- 
luminosity CT AGN up to z^2- 3. Similarly, the planne d 
New X-ray Telescope^ (NeXT; iTakahashi et al.ll2008D . 
scheduled for launch in 2013 will provide imaging and 
spectroscopy in the 5-80 keV energy band with an an- 
gular resolution <1.7' and a spectral resolution of ~ 1.5 
keV. Another focusing hard X-ra y observatory, Simbo l- 
Xi°, is targeted for launch in 2014 (jFerrando et al.ll2004D . 
Simbol-X will perform pointed observations with a field 
of view of ~12' and an angular resolution of ~30". 

Finally, it is important to note that for z~2 the Chan- 
dra and XMM observed energy band of 2-10 keV cor- 
responds to a rest-frame energy of ~6-30 keV, so the 
effects of obscuration are less important. Unfortunately 
even the deepest Chandra data available now only de- 
tect a few photons for the CT AGN candidates at z^2 
(e.g., iTozzi et all l2006l ). thus preventing detailed spec- 

M ore information about EXIST can be found at 

|http :/ /exist. gsfc.nasa.gov/ 

^The NuSTAR website URL is http://www.n ustar.caltech.edu/| 

^ http:/ /www. astro. isas.ac.jp/future/Ne XT/ 

" More details can be found at http:/ /smsc.cnes. fr/SIMBOLX/j 
http;//ixo. gsfc.nasa.gov/ 



tral fitting that could provide a deeper physical under- 
standing of the nature of these sources. The proposed 
International X-ray Observatory-'^^ (IXO) will provide an 
outstanding opportunity to study these highly-obscured 
high-rcdshift sources. As reported by [Alexander et al.l 
12008.) . the IXO will be able to detect thousands of pho- 
tons for the CT AGN detected in the Chandra Deep 
Fields observations for similar, ~1 Msec, exposure times, 
yielding high signal-to-noise spectra for these sources. 
Deep observations at high energies with NuSTAR, EX- 
IST and Simbol-X will provide large samples of heavily 
obscured AGN at 2:~l-3, while the improved sensitivity 
and spectral resolution of the IXO will allow us to study 
in detail the spectra of CT AGN at Z--2-3. 
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